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Abstract. Heteronuclear molecules attract wide attention due to their permanent 
electric dipole moments. Analogous to atoms with magnetic dipoles, the existence of 
nonzero electric dipoles significantly enhance possibilities and mechanisms for control 
and design of quantum degenerate molecule systems with electric (E-) fields. This work 
proposes a vortex creation mechanism inside a condensate of heteronuclear molecules 
through the adiabatic flipping of the axial bias of an analogous E-field Ioffe-Pritchard 
trap (IPT), extending the original protocol due to Ishoshima, et ai, [Phys. Rev. 
A 61, 063610 (2000)] for an atomic spinor condensate inside a magnetic (B-) field 
IPT. We provide both analytic proof and numerical simulations to illustrate the high 
fidelity operation of this vortex pump protocol. We hope our work provides stimulating 
experimental possibilities for the active investigations in quantum degenerate molecule 
systems. 
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1. Introduction 

The impressive progress in atomic quantum gases can be largely attributed to 
technical advances in the control of atomic quantum states utilizing electromagnetic 
interactions. While laser cooling with near resonant light narrows down atomic 
momentum distribution, forced evaporation on trapped atoms favors spatial confinement 
into the quantum degenerate regime. Since the first successful experiment on atomic 
Bose-Einstein condensation, or the controlled production of many bosonic atoms into the 
same quantum ground state of a trapping potential, the field of atomic quantum gases 
has blossomed into one of the leading frontiers in physics, attracting ever increasing 
interest. The ultimate goal along this line of thought is to control the state of an 
arbitrary assembly of an interacting many body system, and either use it to simulate 
other physical systems or to perform quantum computations. 

Significant efforts have been directed toward the above ambitious goal, with the 
immediate pursuit from more recent endeavors targeted at the generation, control, and 
study of cold molecules. These endeavors include important successes such as the BCS to 
BEC cross-over due to tunable interactions induced through B-field magnetic Feshbach 
resonances and the production of quantum degenerate heteronuclear molecules in specific 
low lying ro- vibrational states of their electronic ground states [1] . The latter advance is 
especially illuminating, highlighting the prospects for observing and utilizing anisotropic 
and long range interactions from the permanent electric dipoles. 

Molecule quantum fluids bring in new possibilities due to their richer internal 
degrees of freedom, such as alignment and rotation [2], even after their hyperfine 
spins like those of atomic spinor quantum gases are neglected. In this short note, we 
propose a simple mechanism to generate a vortex state, relying on the internal degrees 
of freedom of a heternuclear molecule, through adiabatically flipping the axial bias of 
an external E- field Ioffe-Pritchard trap (IPT) [3J. The physics of our proposal can be 
analogously mapped onto the successful protocol first proposed [4j and demonstrated 
for magnetically trapped atomic condensates [5j El E] • Our theory hopefully will lead to 
parallel experiments that begin to address new opportunities afforded by the molecular 
internal degrees of freedom. 

This paper is organized as follows. First, we describe our model for dc E- 
field trapping of cold heteronuclear molecules based on the interaction between their 
permanent electric dipoles with a spatially varying static E-field IPT. This is analogous 
to the operation of the famous IPT for magnetic dipoles with an inhomogeneous B-field 
[3]. The heteronuclear diatomic molecule is modeled as a 3D rigid rotor [2]. We then 
discuss the eigen structure of a heteronuclear molecule inside a dc E-field and identify 
states capable of trapped by a local E-field minimum and their associated quantum 
numbers. Next, the rotational properties of such eigen states are studied, which helps 
to illustrate the proposed vortex pump mechanism based on the flipping of the axial 
bias E-field [S El El El El E] • Finally we conclude. 

It is imperative to present this study because this analogy between E-field and 
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B-field and between molecular E-dipole and atomic B-dipole is not a simple one to 
one map. The magnetic dipole of an atom contains both an electronic and a nuclear 
component, respectively, proportional to the electronic spin (S) and the nuclear spin 
(J). Its projection onto the local B-field is approximately proportional to the projection 
of its hyperfine spin F = S + I when the B-field is weak. In contrast, the molecular 
dipole is a fixed constant and pointed along the molecule axis connecting the two nuclei. 
The internal rotational angular momentum of the molecule stays in the perpendicular 
plane of the molecule axis in the absence of an E-field. Despite these differences, our 
study finds that the Isoshima, et al, protocol [U EE [6], originally developed for B- 
field trapped atoms or molecules, remains effective for heteronuclear molecules in a 
E-field IPT. Much of these comparative studies and analogies will become clear from 
the materials presented in the third and fourth sections. 



2. Trapping a heteronuclear molecule with static electric field 

We consider a diatomic molecule composed of two different atomic nuclei and described 
in terms of their center of mass R and relative coordinate r. Neglecting the nuclear and 
electronic spins and their associated interactions, or assume we consider an eigen state 
of the above internal degrees of freedom, the remaining Hamiltonian reads 

p2 p2 

U = ^ + - D(f) ■ E(R) 




2fir 2 dr \ dr J 2/ir 2 (")()> ( ) 

where E(R) denotes the inhomogeneous static E-field, which is slowly varying over 
the molecule size of r, and D(r) is the permanent electric dipole moment operator for 
this specific electronic and internal spin eigen-state [21 [TOl E]. M = M\ + M2 and 
[i = M1M2/ {M\ + M 2 ) are, respectively, the total and the reduced mass of the two 
atoms making up the molecule. 

We adopt the Bohn-Oppenheimer approximation to study the effective trapping of 
the molecule center of mass motion. Within this approximation, the molecular wave 
function for R and r is decomposed into the form ty(r, R) = §(r)ip(R), and we find 



2^dr [ dry 2/xr 2 K ' K ' 



iPm j (?) = V Mj (R)'4>mA?), (2) 



P 



V 2M + Vm ^) = E mMR). (3) 

Now consider for a specific vibrational state, the Eq. ([2]) above then reduces to the 
simple form of a 3D rotator in a dc E-field with a rotational constant B = (h 2 /2/ir 2 ). 
If we further assume the local E-field direction to be the quantization z-axis, then the 
Schrodinger equation for the internal part of the diatomic molecule becomes 

BJ 2 - DE(R) COS0 r l *PmM, 4>r) = V Mj {R)lpMj(0r, <t>r), (4) 
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where we have neglected the dependence of the permanent dipole moment on the 
internuclear distance r. 

Some intuition can be gained for weak E-fields if a simple perturbation theory 
is adopted as is done in Ref. [TT]. We take the unperturbed Hamiltonian to be 
H = BJ 2 with H \J,Mj) = E JMj \J,Mj), where the eigen-energy and eigen-function 
are, respectively, E JMj = BJ(J+1) and | J, Mj), with spherical harmonics | J, Mj) being 
functions of f. The dipole interaction —DE(R) cos8 r is treated as a perturbation. To 
first order the wave function becomes 



^ DE 
if; Mj = \ J,Mj) + 



2 



(j + i) 2 -m; 

B 2(J+ 1)^ (2J + l)(2J + 3) 

DE 1 



' l J + l,Mj) 

T lj TO) 

J 2 - M 2 



~W2j\ (2J-l)(2J + l) |J " 1 ' Mj) ' (5) 
and the corresponding eigen-energy to second order in the external E-field becomes 

V Mj - BJ(J + 1) + - 2(2J _ 1)(2J + 3) - (6) 

For this model system, Mj is a conserved quantity. Assuming the center of mass 
motion is adiabatic with respect to the internal state, the above eigen energy from the 
internal state Vmj then clearly acts as an effective potential due to its dependence on R. 
In particular, we find any state with 1 — 3Mj/J(J+ 1) > is a weak field seeking state, 
i.e., free space E-field traps analogous to the B- field IPT can be constructed. These 
states, of course, are meta-stable, and trapping is possible because of the dynamic 
stability from the Larmor-like precession of the permanent electric dipole along the 
direction of a local E-field. On the other hand, the strong field seeking states with 
1 — 3Mj/J(J + 1) < cannot be used to spatially confine molecules with static fields, 
because it is impossible to construct a local E-field maximum. 

The significant progress gained over the years in trapping and manipulating 
neutral atoms with static B-fields [12] provide important enabling technologies for the 
experimental successes of cold atom physics research. A variety of B-field traps have 
been implemented for various applications [21 E2]- Despite the analogy described earlier 
of the B-field confinement of magnetic dipoles with the E-field confinement of electric 
dipoles (of polar molecules), experimental efforts of trapping and manipulating polar 
molecules with electrostatic E-fields become an active research topic only in recent years 
with the interests on cold atoms broadening into cold molecules. The Meijer's group 
pioneered the research of cooling and trapping of molecules with (3D quadrupole) E- 
fields [13]. They also demonstrated the guiding of cold molecules with a torus shaped 
(2D) E-field hexapole [Hj. With suitable modifications, an analogous IPT of E-field can 
be constructed for polar molecules [31 [121 [HI IE] . 

In Fig. [H we show the calculated dc Stark shift for the Mj = 1 state of heteronuclear 
molecule KRb that is connected to the J = 2 manifold when the E-field is absent. We 
adopt the parameters as from its ground state X 1 J2 v — as measured recently [1] 
with the electric dipole moment D = 1.36 x 10 4 Hz ■ m/V and the rotational constant 
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B = 1.1139 GHz. The calculation is accomplished by a numerical diagonalization of 
the Hamiltonian over a truncated basis of common eigen states for operators J 2 and 
J z . The diagonalization is carried out in the subspace of a conserved Mj [TO], and the 
procedure is found to be rapidly convergent. At the E-field strength we consider, the 
perturbation result ipMj in Eq. (j5J) with Mj = 1 and J = 2 is found to consists of an 
excellent approximation. This particular state as illustrated in Fig. [1] is clearly a weak 
field seeking state. The maximum trap height can become as high as 173 MHz. 

200 I 1 1 1 1 1 1 1 1 1 




Electric Field (kV/cm) 

Figure 1. The Stark shift of the Mj — 1 state for the polar molecule KRb in the 
electronic ground state X 1 v — 0; f° r realistic electrostatic E-trap field strength 

PUCES]- 

3. Rotational transformation of the wave function \J,Mj) 

Before we extend the Isoshima, et al, protocol [4 J from B- field trapped atoms to E-field 
trapped polar molecules, we will briefly review the rotational properties of the spherical 
harmonics \J,Mj). Under a rotation along a unit vector h by an angle 6, the wave 
function changes to 

R{6n)\J,M) = e- m - 3 \J,M). (7) 

Expressed in terms of the Euler angles, the most general 3D rotation takes the following 
form 

R(a, P, 7) = exp(— iaJ z ) exp(—i/3J y ) exp(— ijJ z ). (8) 
According to the representation theory of rotation, we find 
R(aJn)\J,M) = Y,D J M , M (a,l3n)\J,M') 

M' 

= ^ex P HM'«]^, M (/?)exp[-zM 7 ]|J,M'), (9) 
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where 



d J M , M (P) = (J,M'\ex P [-tp.J y ]\J,M) 

= [(J + M)!(J - M)!(J + M')!(J - M')!] 1/2 

x EK" 1 )^ 7 — M' — v)\(J + M — v)\{y + M' — M)!z/!] _1 

V 

2J+M-M'-2v / M'-M+2i' 

(10) 





The value of z/ needs to ensure that the numbers in the factorials stay non-negative. 
For example, when (3 = it and cos | = 0, for d J M , M to be nonzero, we need to let 
2 J + M - M' - 2v = 0. This gives 



(J-M' + z/)! = --(M' + M) 



'J + M-z/)! 



2 

-(M' + M) 



(11) 
(12) 



.2 

In order to assure that both J — M' — u and J + M — v are non-negative, we need to 
make M' + M = 0. As a result we find 

^m'm( 71 ') = (-1) J ~ M <W -Af, (13) 

which gives 

i?(a,7r,7)|J,M) = (-1) J+M exp[iMa] exp[-iM 7 ] | J, -M). (14) 

We want to stress at this point that the angles a, f3, and 7 are parameters used 
to specify arbitrary rotations, and they have nothing to do with the internal state 
angles 9 r and r . To understand the physics behind the vortex generation protocol, 
which is directly related to the geometrical properties of the static E-field, we first 
consider several special cases of unit vectors as rotation axes: for h = X cos + Y sin 0, 
h± = —X sin + Y cos 0, and n q = X sin + Y cos 0, respectively, with X and Y the 
unit vectors along the Cartesian x- and ^/-directions. 6 and denote the polar and 

— * 

azimuthal angles of the center of mass coordinate R. The last case of n q corresponds 
simply to the direction of a IPT E-field in the plane of Z = 0. Using the relationship 
derived above, we find that 

Rn(n)\J,Mj) = exp(-mn- 3)\J,Mj) 

= R,{4> - n/2)Ry(7t)R,(7i/2 - 0)| J, Mj) 

= _ expfiiMrfM-Mj), (15) 

RnM\ J , m j) = exp(-iirn ± • J) I J, Mj) 

= ^((vr/2 + 0) - 7r/2)^(7r)i?,(7r/2 - (vr/2 + 0))| J, Mj) 
= (-l^exp^M^-Mj), (16) 

R Aq (7r)\J, Mj) = exp(-inh q ■ J)\J, Mj) 

= R £ ((n/2 - 0) - 7r/2)^(7r)^(7r/2 - (vr/2 - 0))| J, M,) 
= (_i) Mj exp[— 2iMj0]| J, -Mj). (17) 
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These properties of rotational transformation show that by enforcing a rotation of 
the internal state, the wave function gains an appropriate topologically phase specified 
by the azimuthal angle coordinate of the center of mass (or the molecule) coordinate. 
Provided the internal state flipping of the permanent dipole moment and the center of 
mass motion are both adiabatic, the above results show that different vortical phases 
are generated as in the protocol of Isoshima, et ai, for B-field trapped atomic spinor 
condensates [U [5j [6] . 

4. Vortex creation in a condensate of heteronuclear molecules 

In this section, we will confirm numerically the vortex creation protocol for a condensate 
of heteronuclear molecules in an E-field IPT. Basically, we will simulate the axial E-field 
bias flip and check for the adiabatic conditions of the associated internal state. There 
are two main points we need to watch for. First we need to create the proper vortical 
phase structure on a condensate during the time evolution of the E-field. Secondly we 
must maintain adiabaticity during the flip of the internal state of a polar molecule. 

Considering the first question, assume the initial state corresponds to a polar 
molecule locally aligned along the E-field direction of a E-field IPT, which is essentially 
along the axial z-axis direction close to the z-axis, its local internal state 4>Mj{9r, 4>r) can 
then be described approximately as the eigen state of the system in the local E-field. 
Since Mj is a conserved quantity, this eigen state ipMj(Qr,<f>r) can be expanded by a 
series of | J, Mj) in J, 

^ Mj (0r,<Pr)=J2 C M M j)- ( 18 ) 
J 

The flipping of the E-field bias adiabatically corresponds then to nothing but a rotation 
of the initial state along the unit vector n q (8, <fi) in the transverse plane by an angle of 
7r. According to the rotation properties we discussed before, this flipping of the bias 
then gives 

Rn t (7r)ip Mj (0 r , <p r ) = ^2{—^)' J+Mj Cj exp[—2iMj(j)]\J, -Mj), (19) 

j 

which gains clearly a vortex with a winding number —2Mj. 

Turning to the second question concerning the adiabaticity, we simply can 
propagate the initial wave function ip Mj in real time, simulating the complete process 
of the bias flip. To test the level of adiabaticity and the validity of Eq. (JT9l . we 
consider points (X Q ,Y ,Q) along a circle in the X-Y plane. The absolute value for 
the E-field along the circle is then a constant, although their local directions are all 
different. In fact, for the E-field IPT being considered here, the E-field can be written 
as E = E'(X, —Y, L) with L a constant. Choosing the weak field seeking state ipM J= i 
with J = 2 at E = 3kV/cm Z as considered in Fig. (TJ we simulate the flip protocol over 
various time intervals. The time evolution of the bias E-field is taken to be the same as 
that in the optimal B-field protocol [9], and the E-field in the transverse X-Y plane is 
treated constant. 
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Figure 2. The topological phase gained after numerically flipping the bias E-field, 
that is calculated by taking the inner product of the time evolved state after flip with 
the nominally down polarized (Mj = — 1) state. ■ denotes numerically computed 
data, and the red solid curve represents a linear fit. Both axes are dimensionless in 
units of 2tt. 

In our simulation, we take the E-field at t = as E'L = 3kV/cm and 
E'yJXg + Y§ = 0.3kV/cm. After the flip of the E-bias over W 5 /B (~ 0.1 ms at 
B = 1.1139 GHz), which is a reasonably fast time scale, we calculate the inner product of 
the final state with the eigen-state ipMj=-i (of J = 2) when the E-field E = — 3 kV/cm Z 
is pointed in the downward direction. Using these parameters, the adiabaticity is found 
to be fully satisfied, and the state overlap is essentially 100% except for the extra 
topological phase gained as shown in Fig. [2j A numerical fit gives precisely the slope 
for the phase over being exactly equal to —2, which confirms the high fidelity operation 
of our vortex pump proposal. 

Before concluding, we point out that, like the system of magnetically trapped atomic 
spinors inside a B-field IPT [T71 [9], the quantity of J z — L z is found to commute with 
the Hamiltonian B J 2 — D • E'(X, —Y, L), where E' is the spatial gradient of the E-field 
IPT, and L z is the mechanical angular momentum of the heteronuclear molecule. J is 
the rotational angular momentum of the molecule. If the flip of the E-bias is indeed 
adiabatic, the respective quantum numbers conserve the combination Mj — L z . After 
going through the flip, the internal rotational state is changed from Mj to — Mj, which 
then must be accompanied by an increase of 2Mj to its mechanical angular momentum. 
Thus we see the appearance of a vortex state. 

In conclusion, by analogy with B-field trapping of neutral atoms with magnetic 
dipoles, we study and identify weak field trapping states of polar molecules inside a 
spatially inhomogeneous dc E-field. Further, we suggest that an effective and efficient 
vortex pump protocol can be envisioned for condensates of polar molecules [8l [9], 
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based on the flipping of the axial bias field, originally suggested jlj and experimentally 
demonstrated for atomic spinor condensates inside a B-field IPT El E|- We have 
confirmed that the E-field bias flip remains effective, and the vortex state created 
maintains a vorticity proportional to the Mj quantum number of the trapped molecule 
state. When a diatomic molecule is placed inside a homogeneous E-field, Mj is a good 
quantum number. The weak field trapping states can be associated with very large 
values of Mj, thus the amount of vorticity created could become very significant even 
after a single bias flip. This vorticity possibly could open up a practical approach to 
reach the rapid rotation limit of atomic/molecular quantum gases. 

More generally, we find that this protocol for vortex creation remains effective if 
the diatomic molecule is taken as a symmetric top [18J. Interestingly we find the E- 
bias flip also works for the strong field trapping states provided the polar molecules are 
confined through other means not relying on its permanent electric dipole. Finally, other 
improvements to the B-field bias flip protocol [4] , such as those developed for cyclically 
operated continuous vortex pumping schemes can be analogously extended to the case 
of heteronuclear molecules (3 [9] . 

5. Acknowledgement 

This work is supported by US NSF, NSF of China under Grant 10640420151, and 
NKBRSF of China under Grants 2006CB921206 and 2006AA06Z104. 

References 

[1] Ni K -K, Ospelkaus S, Miranda M H G de, Pe'er A, Neyenhuis B, Zirbel J J, Kotochigova S, 

Julienne P S, Jin D S, Ye J 2008 science 322, 231 
[2] Friedrich Bretislav and Herschbach Dudley 1995 Phys. Rev. Lett. 74, 4623 
[3] Pritchard D E 1983 Phys. Rev. Lett. 51, 1336 

[4] Isoshima T, Nakahara M, Ohmi T, and Machida K 2000 Phys. Rev. A 61, 063610 

[5] Leanhardt A E, Gorlitz A, Chikkatur A P , Kielpinski D, Shin Y, Pritchard D E, and Ketterle W 

2002 Phys. Rev. Lett. 89, 190403 
[6] Leanhardt A E, Shin Y, Kielpinski D, Pritchard D E, and Ketterle W 2003 Phys. Rev. Lett. 90, 

140403 

[7] Isoshima T, Okano M, Yasuda H, Kasa K, Huhtamaki J A, Kumakura M, and Takahashi Y 2007 

Phys. Rev. Lett. 99, 200403 
[8] Mottonen M, Pietila, V, and Virtanen S M M 2007 Phys. Rev. Lett. 99, 250406 
[9] Xu Z F, Zhang P, Raman C, and You L 2008 Phys. Rev. A 78, 043606 
[10] Cohen M, Feldmann Tova, and Kais S 1984 J. Phys. B: At. Mol. Phys. 17, 3535 
[11] Micheli A, Pupillo G, Bitchier H P, and Zoller P 2007 Phys. Rev. A 76, 043604 
[12] Bergeman T, Erez Gidon, and Metcalf Harold J 1987 Phys. Rev. A 35, 1535 
[13] Bethlem Hendrick L, Berden Giel, Crompvoets Floris M H, Jongma Rienk T, van Roij Andre J 

A, and Meijer Gerard 2000 Nature 406, 491 
[14] Crompvoets Floris M H, Bethlem Hendrick L, Jongma Rienk T, and Meijer Gerard 2001 Nature 
411, 174 

[15] Shafer-Ray N E, Milton K A, Furneaux B R, Abraham E R I, and Kalbfleisch G R 2003 Phys. 
Rev. A 67, 045401 



Pumping vortex into a Bose-Einstein condensate of heteronuclear molecules 



10 



[16] van Vcldhovcn Jacqueline, Bethlcm Hendrick L., Schnell Mclanic, and Meijer Gerard 2006 Phys. 
Rev. A 73, 063408 

[17] Zhang P, Jen H H, Sun C P, and You L 2007 Phys. Rev. Lett. 98, 030403 

[18] Herzberg Gerhard 1956 Molecular Spectra and Molecular Structure, I: Spectra of Diatomic 
Molecules (2nd) (D. Van Nostrand Company) 



